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Summary. 3D spectroscopy produces hundreds of spectra from which maps of 
the characteristics of stellar populations (age-metallicity) and internal kinemat- 
ics of galaxies can be derived. We carried on simulations to assess the reliability 
of inversion methods and to define the requirements for future observations. We 
quantify the biases and show that to minimize the errors on the kinematics, age 
and metallicity (in a given observing time) the size of the spatial elements and the 
spectral dispersion should be chosen to obtain an instrumental velocity dispersion 
comparable to the physical dispersion. 
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1 Analyzing the data 

Recently it became possible to derive simultaneously internal kinematics and 
characteristics of the stellar population of galaxies [1]. 

We use a simple parametric procedure [2] to fit the moments of the line- 
of-sight velocity distribution (LOSVD): v, a, and parametrized star for- 
mation history (SFH) containing either single stellar population (SSP) or 
several star bursts, using models computed with PEGASE.HR [3]. Before 
running the inversion procedure one needs to: (a) determine variations of 
the spectrograph's line-spread function along wavelength and over the field 
of the IFU; (b) inject this information into template spectra. We stress that 
the population parameters are constrained by the absorption lines and not 
by the shape of the continuum which may be affected by internal extinction 
or calibration uncertainties. 



2 Validation of the population pixel fitting 

The questions we address are: 

1. Does the method return unbiased estimates of kinematics, age (t) and 
metallicity (Z) at any signal to noise ratio (SNR)? 

2. What are the degeneracies between different parameters? 
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3. What SNR and spectral resolution are required to get a given precision? 

4. Are various models of stellar populations consistent? 

Biases and degeneracy. We have performed extensive Monte-Carlo sim- 
ulations to fit SSP or SFH containing two bursts with SNR ranging from 1 
to 500 pix -1 (at R=10000; pix=0.2A). As long as the grid of models is fine 
enough for performing good interpolation, our method is not biased down to 
SNR=5 pix -1 . The main degeneracy is naturally between age and metallic- 
ity, but Fig 1 presents also the degeneracy between a and Z: a discrepancy 
toward higher Z (sharp lines) is compensated by higher a. This latter de- 
generacy is the strongest coupling between the kinematical and population 
parameters. An error of 1 dex on the metallicity results in an error of about 
25% on the velocity dispersion. This can introduce a significant systematics 
in the determination of the mass-to-light ratio of galaxies. This effect is con- 
siderably reduced if an additive continuum is included in the fit, but in this 
case constraints on the stellar populations are lost. 

Relation between precision and SNR. The precision essentially depends 
on the total SNR integrated over the whole wavelength range (Fig 2): a 
lower SNR per pixel can be balanced by a larger number of pixels. When the 
wavelength range is shortened to exclude the blue region containing line 
and bluer (but keeping the same total SNR) , the precision on kinematics and 
metallicity is not seriously affected, but the precision on age becomes twice 
worse. 
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Fig. 1. Metallicity- velocity dispersion and Age-metallicity degeneracies. Monte- 
Carlo simulations with PEGASE.HR spectra, SNR = 100 pix -1 , a = 50kms _1 



Relation between resolution and precision. The precision on kinematical 
and population parameters depends on (a 1 + of ns ) (Fig 3), where a- lns is the 
instrumental velocity dispersion. The resolution is an important parameter 
for the precise determination of the kinematics. But it has weaker influence 
over the errors on age and metallicity .(See the relations below.) When a < 



Kinematics and Stellar Populations of Galaxies: Validation 



3 



0.4 ai ns , the different minimization strategies that we experimented generally 
become unstable to measure the internal kinematics. 



3 Conclusion 

The method based on the pixel fitting with PEGASE.HR templates can effi- 
ciently constrain kinematics and stellar population. We summarize below the 
relations between the errors on the population characteristics, the total signal 
to noise ratio, and the resolution (for the spectral range 4000 A- 6800A): 

err v « 5 x SNR- 1 x {a 2 + <r? ns ) 2/3 
a x err a ss 5 x SNR- 1 x (a 2 + af ns ) 6 / 5 
err t /t w 54 x SNR- 1 x (a 2 + ct^) 1 / 14 
err z w 12 x SNR' 1 x (cr 2 + er 2 ns ) 1/7 

These relations are useful to select the observational setup: what is the 
optimal compromise between the spectral dispersion and the size of the spa- 
tial elements that minimizes the errors in a given observing time? Consider- 
ing only the sources of uncertainties modeled in our simulations, it appears 
that the best precision on the internal kinematics will be obtained when 
(Ti ns « 0.8(7. For the best precision on the parameters of the stellar popula- 
tion it is preferable to maximize the total SNR with a lower dispersion. Still, 
a single setup with a resolution matching the velocity dispersion is in general 
the good choice. 

The template mismatch due to abundance effects and the uncertainties 
in the modeling of a stellar populations are probably the source of significant 
biases on the parameters of the stellar population. In particular, we can guess 
that if Balmer lines are not in the wavelength range, the determination of 
age will be extremely sensitive to abundance mismatch. We have inverted 
observations of globular clusters [4] and simulated spectra from Bruzual & 
Chariot [5] . For old populations the estimates appear too young and metallic. 
We are investigating the origin of this problem that may be connected with 
the bias found by Prieto et al. [6] in the ELODIE library. 
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SNR [per -pixel, for R-10000, pixel-1 km/s, range=400-680nm] 
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Fig. 2. Error on population and kinematical parameters as a function of the total 
SNR. The upper x axis shows SNR per pixel for A range from 400 to 680 nm and 
dispersion of 0.2A/pix. The curves noted 'BC03' are representing simulations with 
Bruzual & Chariot (2003) models, and those noted 'Rose GC are observation of 
globular clusters by Rose et al. (2005), all analyzed with PEGASE.HR templates. 
Errors are proportional to (total SNR) -1 , see the formula in the text. 
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Fig. 3. Errors on the population and kinematical parameters as a function of the 
resolution. The instrumental dispersion ranges from <Ji n3 — 12 km/s to Oi ns = 
250 km/s (R=500 to 10000). Errors are increasing with the effective broadening 
yj (tj 2 + of ns ) according to the formula given in the text. 



